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Abstract. Numerical hydrodynamical modelling of supernova- driven shell formation is performed with a purpose 
to reproduce a giant HI ring (diameter 1.7 kpc) in the dwarf irregular galaxy Holmberg I (Ho I). We find that 
the contrast in HI surface density between the central HI depression and the ring is sensitive to the shape of the 
gravitational potential. This circumstance can be used to constrain the total mass (including the dark matter 
halo) of nearly face-on dwarf irregulars. We consider two models of Ho I, which differ by an assumed mass of 
the dark matter halo Mh- The contrast in HI surface density between the central HI depression and the ring, as 
well as the lack of gas expansion in the central hole, are better reproduced by the model with a massive halo of 
Mh = 6.0 X 10^ Mq than by that with a small halo of Mh = 4.0 x 10® M©, implying that Ho I is halo-dominated. 
Assuming the halo mass of 6.0 X 10^ Mq, we determine the mechanical energy required to form the observed ring 
equal to (3.0 ± 0.5) x lO''^ ergs, equivalent 300 ± 50 Type II supernovae. The inclination of Ho I is constrained to 
15° — 20° by comparing the modelled HI spectrum and channel maps with those observed. 
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^ 1. Introduction produce a ring-like morphology of the distribution of the 

atomic hydrogen in Ho I. We show that modelling of the 

^ The interstellar medium of gas-rich dwarf irregular galax- HI spectrum, HI channel maps, and the radial HI profile 

■ ies (dirr) is usually dominated by features described as may constrain the inclination angle and consequently the 

shells, holes, or rings. The number of such structures varies dark matter content of Ho I. 

from a few (IC 10, Wilcot s fc M iller [1998|) to several tens in Sect. El the origin of a ring-hke HI morphology in 

(Holmberg II, Puche et al. [199211 and their sizes are within Ho I is discussed. In Sect.Elthe numerical hydrodynamical 

several hundred parsecs. However, there are a few dirr's model for simulating the supernova-driven shell dynamics 

the HI morphology of which is totally dominated by a sin- jg formulated. The results of simulations are presented in 

gle ring str ucture of the size comparable to their optical Sect. H and Sect. Our conclusions are summarized in 

extent (Ott 1199"^ . Among them is Holmberg I, a mem- Sect. El 
ber of the M81 group of galaxies (distance ~ 3.6 Mpc). 
Ott et al. H200ip have found that most of its HI content 

(75%) is localized within a giant ring of 1.7 kpc diameter. 2. Origin of the peculiar HI morphology in Ho I 

They suggested that strong stellar winds and supernova Fig. [H shows the integrated HI emission of Ho I with the 
(SN) explosions might be responsible for the peculiar HI kinem atical major and minor axes overlaid (Ott et al. 
morphology in Ho I. IMTTj) . Most of the HI content is obviously concentrated 
In the present paper we test this hypothesis via de- in a ring, the morphological center of which is offset by 
tailed numerical hydrodynamical modelling. Two model 0.7 kpc from the dynamical center of the galaxy. The con- 
galaxies are considered: a fast rotating galaxy with a mas- trast in column density between the central HI depression 
sive halo and a slowly rotating one with a small halo. We and the ring is about a factor of 15. In principle, kpc- 
attempt to determine the amount of SNe required to re- sized ring structures in the gas component of disk galaxies 
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Fig. 1. Integrated HI emission of Ho I. In the bottom left 
comer we show the half-power beam (8'.'2 x 7.0'.'). The 
grayscale represents the HI column density. The data have 
been obtained from the Very Large Array observations 
(Ott et al. nUUT^ . 



can be caused by a number of phenomena. Among them 
are resonant rings ('Buta lTUS^ . coUisionally induced rings 
(Lynds & Toomre I1976|l , the infall of massive gas clouds 
(Tenorio-Tagle et al. I1987|) . gamma ray bursts (Efremov 
et al. I1998|l . or multiple supernova explosions. In case of 
Ho I, the first two scenarios can be ruled out because of the 
absence of spiral structure, the apparent absence of lower- 
mass companions along the projected minor axis of Ho I, 
and the lack of expansion motions in the ring. The latter 
three scenarios are plausible. However, the HI ring in Ho I 
shows evidence for a supernova-driven origin, since the B- 
band image of Ho I indicates that there is a pronounced 
concentration of young blue stars within the ring (Ott et 
al. 12001(1 . and the most massive stars of this population 
might have been a plausible source of energy input over 
the past 50 Myr in our model. In the following sections the 
supernova-driven mechanism is numerically investigated. 
Other possibilities of HI ring formation in Ho I will be 
discussed in a subsequent paper. 

3. Numerical model 

3.1. Initial conditions 

We consider an axisymmetric model galaxy consisting of 
a rotating gas disk, a rigid stellar disk, and a rigid, spher- 
ically symmetric halo. Extensive radio and optical obser- 
vations of Ho I by Ott et al. H2()()l|l provided important 
observational constraints on the total mass and density 
distribution of different components of our model galaxy. 
The HI mass of Ho I is 1.1 x 10^ Mq, with the total gas 
mass (A/hi + A/hc + Mu^) amounting to 1.5 x 10^ Mq. 
Using a B-haud luminosity Lb = 1-0 x 10® Lq and a mass- 



to-light ratio Ms/Lb — 1, Ott et al. (2001) derived an es- 
timate of the luminous stellar mass, Ms = 1.0 x 10® Mq, 
in Ho I. For computational purposes, the density profile 
of the stellar component is chosen as: 



Ps = Ps sech(z/zs) exp(-r/rs), 



(1) 



where p^ is the stellar density in the center of the galaxy, 
and Zs and rg are the vertical scale height and radial scale 
length of the stellar component, respectively. The radial 
scale length Tg can be estimated from the Jc-band radial 
surface brightness profile of Ho I (Ott et al. ^001). If the 
mass follows the light at the longer wavelengths, it results 
in Tg ~ 1.5 kpc. We have adopted a value of 300 pc for the 
vertical scale height Zs, which is typical for dwarf irregu- 
lars. With and Zg being fixed, the stellar density at 
(z, r) = is varied so as to obtain the measured luminous 
stellar mass of Ms = 1.0 x 10® Mq within the computa- 
tional domain. This results in p^ « 0.02 Mq pc^'^. 

The dark matter content of Ho I is uncertain. We 
assume that the density profile of the halo can be ap- 
proximated by a modified isothermal sphere (Binney & 
TremaineUHHIll 



Ph = 



Pw 



1 + (r/rh)2 



(2) 



where the central density p^o and the characteristic scale 
length Th were given by Mac Low & Ferrara p999|l and 
Silich & Tenorio-Tagle lPM|l : 



Pho = 6.3 x 10^ 



Th = 0.89 X 10" 



Mq 

Mh 

Mq 



-1/3 



h-^^^ Mq kpc" 



1/2 



(3) 



(4) 



Here, h 
100 km s 



is the Hubble constant in units of 
Mpc^^ and Mh is the total halo mass. 
We adopt h = 0.65 throughout the paper. 

We consider two different initial models that attempt 
to represent the appearance of Ho I before the ring has 
been formed. 

a) Fast rotating galaxy with a massive halo (hereafter 
model 1). In model 1 the gravity force of a massive halo 
of 6.0 X 10^ Mq and a stellar disk of Mg = 1.0 x 10® Mq 
is balanced by a fast rotating gas disk. We note that the 
adopted halo mass of 6.0 x 10^ Mq is in fact the total halo 
mass of Ho I. Particularly, the halo mass confined within 
the HI diameter of Ho I (5.8 kpc) is 3.3 x 10® Mq, which 
agrees with the corresponding upper limit of 3.1 x 10® Mq 
derived by Ott et al. We vary the rotation curve 

until the initial gas surface density distribution resembles 
that of Ho I at radii that are unperturbed by shell expan- 
sion, r > 2 kpc. We also make sure that the total HI mass 
within the computational area is equal to the measured 
one (A/hi — 1-1 x 10® Mq). For the adopted halo mass, 
the central density and characteristic scale length become 
Pho — 0.04 Mq pc"'^ and — 0.55 kpc, respectively. The 
adopted halo density profile and the gas rotation curve are 
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Fig. 2. The radial density profiles of the halo and the ini- 
tial rotation curves for model 1 (filled triangles and filled 
squares, respectively) and model 2 (open triangles and 
open squares, respectively). 
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Fig. 3. The equilibrium radial surface density profiles and 
the vertical scale height of the g function of galac- 

tocentric radius for model 1 (filled squares and filled tri- 
angles, respectively) and model 2 (open squares and open 
triangles, respectively). 



shown in Fig. [21 by the filled triangles and filled squares, 
respectively. 

b) Slowly rotating galaxy with a small halo (hereafter 
model 2). In model 2 the gravity force of a small halo of 
Mh = 4 X 10^ Mq and a stellar disk of = 1.0 x 10^ Mq 
are balanced by a slowly rotating gas disk. For the adopted 
halo mass, the central density and characteristic scale 
length become pho = 0.1 Mq pc~^ and rh = 0.14 kpc, 
respectively. Thus, in a slowly rotating galaxy the halo 
density has a rather cuspy profile as shown by the open 
triangles in Fig. |21 



3.2. Equilibrium solution 

Once the stellar and halo density profiles are fixed, we 
can proceed to obtain the initial gas density distribution 
by solving the steady-state momentum equation, 



-\/P = -(v ■ V)v 
P 



(5) 



where P is the gas pressure, p is the gas density, $h and 
$s are the gravitational potentials of the halo and stellar 
disk, respectively. We assume the gas distribution to be 
initially symmetric with respect to the rotation axis and 
equatorial plane, which naturally implies the use of cylin- 
drical coordinates {z,r) for solving Eq. ||SJ). If the gas is 
initially isothermal, then Eq. (jsj can be transformed to 
the following set of equations: 



d Inp 

dr 
d Inp 

dz 



1 V. 



1 



1 

dr ' 



1 /^^ \ 1 rf$s 
(7^\ ^Jz CT^ dz 



(6) 
(7) 



Here, a is the gas velocity dispersion. Ott et al. (2001) 
found an average value of cr = 9 km s~^ for Ho I. The 



halo gravity force can be expressed as: 
V$h = V$h(r,) = 'inCpwrl/rl x 



[r*/rh — arctan(r,/rh)] 



(8) 



where r: 



' + z^ is the radial distance from the galactic 



center, r is the galactocentric radius, z is the height above 
the midplane, and ej.^ = r,/r* is a unit vector. The r- and 
z-components of the halo gravity force (V$h)r, (V'i>h)z, 
can be found by projecting V$h on the corresponding 
coordinate planes. The stellar gravitational potential is 
found by numerically solving the Poisson equation: 



= AnGps. 



(9) 



Eqs. © and O are finite-differenced on a mesh of 
350 X 300 grid points in radial and vertical directions, re- 
spectively, and solved using the lower-upper decomposi- 
tion method to obtain the equilibrium gas density distri- 
bution in a computational domain of 3.5 kpcx3.0 kpc. The 
equilibrium radial surface density profiles and the expo- 
nential scale heights of the gas for both models are shown 
in Fig. 13 In both models the gas has an exponential ra- 
dial profile, which is in agreement with observations of the 
HI distribution in many dirr's (Taylor et al. 11994(1 . There 
is also a pronounced difference between the models with 
massive and small halos, the latter has a thicker gas disk. 
Indeed, in model 2 the cuspy halo mainly pulls the mat- 
ter toward the nucleus, while in model 1 the halo density 
has a somewhat shallower distribution, thus contributing 
significantly to the vertical gravity pull at larger radii and 
subsequently assisting to maintain a thinner equilibrium 
gas configuration. 

3.3. Computational techniques 

We use the ZEUS-2D numerical hydrodynamical code in- 
corporating a rotating gas disk initially in equilibrium in 
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the fixed stellar and halo potentials. A usual set of hydro- 
dynamical equations in cylindrical coordinates is solved 
using the method of finite-differences with a time-explicit, 
operator split solution procedure described in detail in 
Stone & Norman 1)1992(1 . We have implemented the cooling 
curve of Spaans & Norman (1997) and Wada & Norman 
(2001) for a metallicity of one tenth of solar, which is close 
to the value of Z = ^0/12 derived by Miller & Hodge 
p99(i|l for Ho I. The use of cooling function simplifies the 
implementation of cooling by collecting the effects of vari- 
ous coolants. The cooling processes taken into account are: 
(1) recombination of H, He, C, O, N, Si, and Fe; (2) col- 
lisional excitation of HI, CI-IV, and OI-IV; (3) hydrogen 
and helium bremsstrahlung; (4) vibrational and rotational 
excitation of H2 ; (5) atomic and molecular cooling due to 
fine-structure emission of C, C-I-, and O and rotational 
line emission of CO and H2. We note that using an equi- 
librium cooling function, oversimplified particularly in the 
intermediate-temperature range of T < 10"* K, is a com- 
mon practice in modelling the SN-driven shell dynamics 
(see e.g. Mac Low & Ferrara 119991 Murakami & Babul 
I1999|l . Such an oversimplification is justified by the fact 
that most amount of the thermal energy of the shell is lost 
in the high-temperature [T ^ 10^ — 10^ K) layers, while 
the energy lost at the low-temperature end contributes 
less to the total energy balance. Variations of the cool- 
ing rate connected with the fractional ionization are im- 
portant mostly to determine the cooling time of the low- 
and intermediate-temperature layers. However, because at 
these temperatures and the corresponding densities the 
cooling time is very short in comparison to the dynamical 
time, the overall dynamics is insensitive to details of the 
cooling rate at T < 10^ K. 

We use an empirical heating function tuned to bal- 
ance the cooling in the background atmosphere so that 
it maintains the gas in hydrostatic and thermal equilib- 
rium and may be thought of as a crude model for the 
stellar energy input. The disadvantage of this method is 
that the time-independent heating becomes strongly un- 
balanced in the hot regions filled with supernova ejecta. 
This may result in a spurious energy input at later phases 
of the shell expansion when it occupies a large computa- 
tional volume. To minimize this numerical effect, we have 
introduced a tracer field function defined by Yabe & Xiao 
H1993|l . The tracer field is advected with the same algo- 
rithm as the gas density in order to follow the hot gas 
ejected by the supernovae. Heating is prohibited in the 
regions where the tracer field is present, thus largely re- 
ducing the effects of spurious heating of the bubble's in- 
terior by the time-independent heating function. This is 
physically justified since most of the heating in the warm 
interstellar medium comes from the photoelectric heat- 
ing of polycyclic aromatic hydrocarbon molecules (PAHs) 
and small grains, which will be either evaporated or highly 
ionized in the hot (> 10^ K) bubble filled with supernova 
ejecta. Cooling and heating are treated numerically using 
Newton-Raphson iterations, supplemented by a bisection 
algorithm for occasional zones where the Newton-Raphson 



method does not converge. In order to monitor accuracy, 
the total change in the internal energy density in one time 
step is kept below 15%. If this condition is not met, the 
time step is reduced and a solution is again sought. 

The energy of supernova explosions is released in the 
form of thermal energy in the central region with a radius 
of four zones. We use a constant wind approximation de- 
scribed in detail in Mac Low & Ferrara l)1999|l . There is 
a certain degree of freedom in how to distribute 10^"'^ ergs 
of mechanical energy released by a single supernova. We 
decided to convert it totally into a thermal energy, since in 
the present simulations we deal with large stellar clusters 
with hundreds of supernovae. With such an amount of SN 
explosions, the surrounding ISM will be quickly heated 
and diluted, making radiative cooling inefficient. We ad- 
mit that at the initial stages of the stellar cluster evolution 
part of the energy of SNe may be radiated away due to 
the radiative cooling. Hence, our numerical simulations 
provide a lower bound on the number of SN explosions 
needed to create the observed ring. We choose the energy 
input phase to last for 30 Myr, which roughly corresponds 
to a difference in the lifetimes of the most and least mas- 
sive stars capable of producing SNe in a cluster of simul- 
taneously born stars. 

4. Multiple Supernova Explosions 

4.1. Fast rotating galaxy with a massive halo 
(model 1 ). 

We vary the energy input of the starburst in terms of SN 
explosions in order to reproduce a ring-like gas morphol- 
ogy of Ho I. We start by showing in Fig. ^ the temporal 
evolution of the distribution of the gas volume density for 
an energy input equivalent to 300 successive SN explo- 
sions. SNe generate a supersonically expanding wind that 
compresses the surrounding gas, thus creating a bubble 
filled with hot ejected gas. At the end of the energy in- 
put phase (t = 27 Myr) the shell loses its spherical form 
transforming into a prolate spheroid. As a result, tangen- 
tial motion vortices develop in the hot ejected gas, which 
in turn trigger Kelvin- Helmholtz instabilities in the com- 
pressed expanding gas layer of the shell. This results in 
a characteristic ripple-like form of the shell. At approxi- 
mately 40 Myr the shell breaks out of the disk, Rayleigh- 
Taylor instability ensues from the shell acceleration, cre- 
ating a characteristic bubble-and-spikc morphology, first 
mentioned in numerical experiments by Mac Low et al. 
(1989J. Soon after the shell stalls and starts to collapse 
in radial direction, simultaneously pushing the hot ejecta 
to a higher altitude. In fact, part of the galaxy ISM will 
eventually be lost by Ho I, because its velocity exceeds 
the local escape speed. We compute the fraction of the 
total gas mass of Ho I lost by such an outfiow at the time 
when the ring has collapsed [t — 72 Myr) by summing 
up the mass of the gas expanding with velocities exceed- 
ing those of the local escape speed. The fraction is low, 
usually within 1 — 2%. At t ^ 54 Myr the shell appears 
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Fig. 4. Temporal evolution of the gas volume density distribution in model 1, with the energy input equivalent to 300 
successive SN explosions. The velocity field is normalized by the local escape speed. In general, the escape speed is 
a function of both the galactocentric distance and the height above the midplane. Normalizing the velocity field by 
the value of local escape speed helps to differentiate the gas that can potentially escape the galaxy. The contour line 
delimits the region filled with SN-ejected material. The grey-scale bar is in log units of g cm~^. 



as a ring if viewed face-on, and its diameter matches best 
that of HI ring in Ho I. At later times of the evolution 
{t > 70 Myr), the central depression is totally filled and 
Ho I appears as a dwarf irregular, with a declining average 
gas surface density profile. Simulations indicate that for a 
smaller starburst capable of producing only 100 - 150 SNe, 
the shell never breaks out of the disk and hardly expands 
to 1.7 kpc diameter, indicating that 100-150 SNe are not 
sufficient to account for the size of HI ring in Ho I. For a 
larger starburst of 450-600 SN, the shell breaks out after 
approximately 35 Myr. However, the radial expansion of 
the shell proceeds to a larger radius than that of Ho I's 
ring, mainly because of the higher initial momentum ac- 
quired by the shell. 



The total Ha luminosity of Ho I is 4.3 x 10^^ ergs s~^ 
(Miller & Hodge 1996), which implies a star formation 
rate of only 0.004 Mq yr"-'^. We searched for dense gas 
clumps in the expanding shell that could be Jeans unsta- 
ble. We assume that the gas in a restricted computational 
domain becomes Jeans unstable if tg < Ar/vs, where tg 
is the free-fall time, Ar is the size of a computational cell, 
and Vs is the local speed of sound. We also combined the 
neighbouring cells to search for a larger scale (< 40 pc) 
Jeans unstable gas clumps. Our simulations show that on 
the scales of a giant molecular cloud 10 — 40 pc) the 
shell is stable against self-gravity. We note however that 
our simulations were restricted in resolution to 10 pc, and 
thus we did not follow possible gravitational instability 
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Fig. 5. Model HI spectrum of Ho I obtained in model 1 
for three assumed angles of inclination, i = 12°, 40°, and 
70°. The measured HI spectrum is plotted with the filled 
squares. 

on smaller scales (< 1 pc) corresponding to the subunits 
and the cores of molecular clouds that can be in princi- 
ple Jeans unstable, while the clouds as a whole are stable 
(Scoville & Sanders nWjl . 

Furthermore, we attempt to reproduce the HI spec- 
trum of Ho I, shown in Fig. [S] by the filled squares (Ott 
et al. I2001|l and scaled to Ho I's systemic velocity of 
141.5 km s~^. Its shape is remarkably well described by a 
Gaussian with a FWHM of 27.1 km s~^. In order to model 
the HI spectrum, we assume that the gas is thermalized 
in each computational cell, with the velocity dispersion 
controlled by the local gas temperature. We also assume 
that the HI distribution in Ho I is optically thin, which 
is justified considering its rather low HI surface density, 
5^Hi < 10 Mq pc^^. We use the following conversion for- 
mula that links the HI mass per velocity channel Mhi/Av 
with the HI flux density Sm (Binney & Merrifield [T^H|) : 



^ = 2.35 X 10^ 
Av \Mpc 



km t 



(10) 



where D is the distance to Ho I. 

We have constructed the model HI spectrum of Ho I 
for different inclination angles. We can surprisingly well 
reproduce the measured HI spectrum for a small inclina- 
tion angle, i = 12°, represented by the solid line in Fig.O 
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Fig. 6. The relative discrepancy (open circles) between 
the model and measured HI spectrum of Ho I as a func- 
tion of assumed inclination obtained for model 1. The 
flux-weighted relative discrepancy (filled squares) is also 
provided to de-emphasize the low-mass wings of the HI 
spectrum. 

the spectrum for smaller inclinations practically coincide 
with the i = 12° one. The model HI spectrum is taken at 
i = 54 Myr when the face-on (or zero-inclination) gas dis- 
tribution best reproduces that of Ho I. A small difference 
between the peak positions in the measured and model 
profiles most probably indicates that there is a small off- 
set between the morphological center of the ring and the 
dynamical center of Ho I, a feature supported by the mea- 
surements of its HI velocity field (Ott et al. l2U(Jl|l . Another 
possibility is that the starburst responsible for creating the 
HI ring is not exactly located in the midplane, but instead 
closer to the near side of the galaxy. 

As we employ progressively higher inclination angles, 
the agreement between the model and the measured HI 
spectrum deteriorates. In Fig. we plot the relative dis- 
crepancy between the model and measured HI spectrum 
averaged over 70 velocity channels, starting at —35 km s""'^ 
and ending at -|-35 km s^^. The flux- weighted relative 
discrepancy is also calculated to de-emphasize the low- 
mass wings of HI spectrum. The best agreement is obvi- 
ously found for low inclination angles, i ^20°. The near- 
Gaussian shape of the HI spectrum indicates that the 
thermal+turbulent motions dominate the line-of-sight ro- 
tation+expansion motions in Ho I. 

As is seen from Fig. 1^1 modelling of the HI spectrum 
provides an upper limit to the inclination of Ho I. A lower 
limit can be obtained by constructing the model HI chan- 
nel maps and comparing them with the observed ones pub- 
lished by Ott et al. H2001|l . The latter show an interesting 
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tendency: the complete ring appears only in a few channel 
maps near the systemic velocity, while the other channel 
maps show only part of it. Specifically, the part of the ring 
that is faint on the approaching side seems to be luminous 
on the receding side of Ho I and vice versa. We smoothed 
the [128-131] kms-i channel map of Ho I (see Fig. 3 of Ott 
et al. l2001,l to a resolution of 20'.'6 x 21'.'2 and found that 
the relative amplitude of the azimuthal variations in the 
HI intensity along the ring is a factor of 2.5. Furthermore, 
we constructed model channel maps at t = 54 Myr for 
inclinations of 0° - 20°. We find that the observed varia- 
tions in HI intensity along the ring can be obtained only 
if the inclination is constrained to the range 15° - 20°. As 
an example, we plot the model channel maps constructed 
at t = 54 Myr for an inclination of 15° in Fig. [7| It is 
readily seen that the area near the central position shows 
virtually no emission throughout all channels, a feature 
also reported by Ott et al. H2001|l . The asymmetry of the 
model channel maps indicates that the gas is participat- 
ing in both, rotation and expansion motions, the latter 
dominating at higher altitudes (see Fig. 0J). The relative 
amplitude of azimuthal variations in the maximum HI in- 
tensity along the ring in the [-13.5; -10.5] km s^^ channel 
map, which is equivalent to the observed [128-131] km s^^ 
channel map of Ott et al. ()^()01|l . is a factor of 2. This 
value decreases, as we take progressively lower inclination 
angles and vanishes at zero inclination. Hence, if Ho I is 
halo dominated as assumed in model 1, then its inclina- 
tion angle should be limited to a range of 15° - 20°. The 
lower right panel in Fig.Qshows the modelled, integrated 
HI emission. Visual inspection shows that the size of the 
ring and the HI column density distribution are similar to 
those in Fig. indicating that the HI ring can indeed be 
a result of multiple SN explosions. 

The rotation curve of model 1 (the filled squares. 
Fig. EJ agrees well with the observed HI rotation curve 
of Ho I for i = 15° (see Ott et al. OnHTl Fig. 9) only 
in the outer regions at r > 1.2 kpc, which are unper- 
turbed by the shell expansion. In the inner perturbed re- 
gions the observed velocities exceed those of model 1 by 
5-10 km s^^. We note here that this discrepancy might be 
due to an offset location of the dynamical center of Ho I 
with respect to the geometrical center of the HI hole (see 
Fig-Di the effect that can hardly be modelled in our two- 
dimensional axisymmetric simulations. As a result, the 
mass redistribution due to the off-center SN explosions 
might have brought a considerable fraction of the fast ro- 
tating gas towards the galactic center. Indeed, adopting 
the observed rotation curve of Ho I as the initial rotation 
curve of model 1 would result in the equilibrium gas distri- 
bution that totally lacks any gas in the central 1 kpc radius 
irrespective of the halo mass we use; a realistic equilibrium 
gas distribution cannot be obtained in this case, implying 
that the observed rotation curve of Ho I is strongly per- 
turbed in the inner regions. 

Finally, in Fig. |Slwe present the azimuthally averaged 
radial gas distribution at three different phases of the shell 
expansion, taking the morphological center of the pro- 




Fig. 7. Model channel maps of Ho I at the assumed in- 
clination of i = 15°. The lower right panel shows the in- 
tegrated HI emission obtained in model 1. The grey-scale 
wedge is in units of cm~^ and the velocity channels are in 
units of km s~^. 
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Fig. 8. Azimuthally averaged radial gas distribution ob- 
tained with model 1 at three different phases of the shell 
expansion. The measured radial gas distribution of Ho I 
is plotted with the filled squares where a possible contri- 
bution of He and H2 to the total gas mass is taken into 
account (phi+Ho+Hj/phi = 1.4). 

jected shell as the origin. An inclination i = 15° is as- 
sumed. The best agreement with observations is found at 
t = 55 ± 5 Myr, which implies that the shell in Ho I has 
attained its largest size at the present time. The lack of ex- 
pansion motions reported by Ott et al. H2()()l§ is thus not 
due to the low inclination of Ho I, but rather an intrin- 
sic property. The measured contrast ('^ 15) in gas surface 
density between the central depression and the ring is well 
reproduced. We varied the total number of SN explosions 
(A'sn) in order to see how the energy input can influence 
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our conclusions. As mentioned above, for iVsN < 150 the 
shell never expands to its present diameter of 1.7 kpc, 
thus ruling out a small starburst as the origin of the HI 
ring in Ho I. On the other hand, for TVgN > 450 the pro- 
jected shell appears too sharp-peaked at the time when 
it has reached a 1.7 kpc diameter, very similar to the ra- 
dial gas profile for A^sn = 300 at 27 Myr (see Fig. |H1 
dashed line) . Note that a possible depletion of HI column 
density (Ahi) due to the ionization by an external UV 
background with ^ 10^ Lya photons cm~^ s~^ Hz~^ ac- 
counts for AA^Hi 3 X 10^^ N^^ cm^^ only, which is 
less than 14% at t > 54 Myr. Hence, we find that the 
energy output equivalent to 300 ± 50 SNe can be respon- 
sible for producing the HI ring-like distribution of Ho I. 
If we assume a star formation efhciency of 5% — 10% and 
a Salpeter IMF, the mass of a parent gas cloud must be 
arotind 1.5 — 3.0 x 10^ Mq, which is not extreme. 

4.2. Slowly rotating galaxy with a small halo 
(model 2). 

Our numerical hydrodynamical simulations reveal consid- 
erable differences in the dynamics of the SN-driven shell 
for the case of a slowly rotating galaxy with a small, cuspy 
halo (model 2), as compared to the fast rotating galaxy 
with a massive, shallow halo (model 1). First, we find 
that the energy input needed to create the observed ring- 
like morphology is a factor of 2 lower in model 2 than in 
model 1. The radial gas surface density profiles in both 
models are similar; however, the exponential scale height 
of the gas is higher in model 2 (see Fig. which results 
in a smaller gas volume density. This fact, along with a 
much lower inward gravity pull by a small halo, renders 
the isotropic expansion of the shell much more easy, but 
on the other hand impedes break through. Fig. ^ shows 
the time evolution of the distribution of the gas volume 
density for an energy deposition equivalent to 150 SN ex- 
plosions. The shell evidently never breaks out of the disk. 
This is mainly due to the thicker vertical gas distribution 
and the lower energy input in model 2, as compared to 
model 1. On the contrary, the lack of gas expansion in the 
central HI depression in Ho I reported by Ott et al. 120011 
implies that its shell has broken out of the disk. 

Applying the same analysis as in Sect. 14. ll wc find the 
best agreement with the HI spectrum and channel maps of 
Ho I for inclination angles i = 25° - 30°, in contrast to i ~ 
15° - 20° for model 1. Thus, numerical models with differ- 
ent halo masses predict different angles of inclination for 
Ho I. Comparison with the inclination angle derived from 
observations could therefore help to constrain the halo 
mass of Ho I. Unfortunately, Ott et al. (I|2001il were not able 
to get stable results for the inclination. This implies that 
Ho I has a rather low inclination, which favours model 1 
as the prime choice of Ho I. Furthermore, in case of Ho I 
we do not see any HI along the line of sight to the cen- 
tral hole down to a column density of ^ 6.0 x 10^^ cm^^, 
or a surface density of 0.45 Mq pc~^, a feature which 
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Fig. 10. Azimuthally averaged radial gas distribution ob- 
tained with model 2 at three different phases of the shell 
expansion. The measured radial gas distribution of Ho I is 
plotted with the filled squares. 

is much better reproduced in model 1 than in model 2. 
For instance, the gas surface density near the central de- 
pression in Fig. 1101 never falls below 6 Mq pc~^. Indeed, 
the shell has not yet broken out of the disk in model 2, 
thus contributing substantially to the line-of-sight gas col- 
umn density near the central depression. This tendency is 
clearly seen in Fig. 1101 which shows the azimuthally av- 
eraged radial gas distribution at three diff'erent phases of 
the shell expansion, taking the morphological center of the 
projected shell as the origin. An inclination of i = 25° is 
assumed. As is seen, the surface density of the gas near 
the center always exceeds that of Ho I shown by the filled 
squares. 

5. Discussion 

Ho I has a ring-like HI distribution, which is moderately 
non-axisymmctric with respect to its morphological cen- 
ter (see Fig.^. The ellipticity of the ring makes it rather 
difficult to draw any quantitative comparison with the ob- 
servations on the basis of two-dimensional axisymmetric 
simulations. However, the basic qualitative features of the 
gas distribution and velocity structure in Ho I, such as 
the size of the ring, the contrast in the gas surface den- 
sity (hereafter, Sg) between the central depression and the 
ring, and the lack of gas expansion in the central hole can 
be modelled in the two-dimensional case. 

The lack of gas expansion in the central depression 
reported by Ott et al. f(2001| argues that the shell has 
suffered a blowout. Otherwise, one would expect the shell 
to be still in the Sedov expansion phase with some pilcd- 
up gas in the z-direction. However, we do not see any HI 
along the line of sight to the central hole down to a col- 
umn density of ~ 6.0 x 10^^ cm~^. The measured contrast 
in Sg between the central depression and the ring (a fac- 
tor of 15) is also indicative of a blowout, since the 'limb 
brightening' effect may account for a factor of 2 difference 
in Eg between the projected rim of a three-dimensional 
shell and its center. Hence, the occurrence of a blowout 
(implied by the measured contrast in Eg between the cen- 
tral depression and the ring) at the time when the shell 
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Fig. 9. Time evolution of the gas volume density obtained with model 2, with an energy input equivalent to 150 
successive SN explosions. The velocity field is normalized to the local escape speed. The contour line delimits the 
region filled with SN-ejected material. The grey-scale wedge is in units of g cm^"^. 



has expanded out to its present size of 1.7 kpc is the main 
criterion for the selection between the two models. 

As discussed in Sect. l4.ll and l4.2l the blowout scenario 
of Ho I is better reproduced in model 1 at i = 55 ± 5 Myr. 
At a first glance, the radial gas distribution in model 2 at 
73 Myr seems by eye to better fit the overall observational 
data. However, it is not only the measured size of the ring 
(1.7 kpc), but also the measured contrast in Sg between 
the central depression and the ring (a factor of 15) that 
are not reproduced for this case. We would like to note 
here that the estimated decrease in HI column density 
near the central depression caused by the external UV 
background is approximately 2% at i > 45 Myr in model 2. 
Hence, the external UV radiation field cannot reconcile the 
measured contrast of 15. Although the overall radial gas 
distribution in model 2 at 73 Myr fits better the measured 
profile, it fails to reproduce both the size of the ring and 
the measured contrast in the gas surface density between 
the central depression and the ring. 



On the other hand, our best model is not optimized 
to fit the observations. At the time when the ring has at- 
tained its present size of 1.7 kpc, the maximum gas surface 
density exceeds that of Ho I by approximately 5 Mq pc~^ 
(see FiglHl the solid line). At the same time, the width 
of the modelled ring is smaller than that of the observed 
ring. This is most probably due to the two-dimensional 
axisymmetric nature of our simulations attempting to re- 
produce a moderately non-axisymmetric structure, or az- 
imuthal averaging of the inclined ring with vertical walls. 
Three-dimensional non-axisymmetric modelling of Ho I is 
required for a fine tuning of the model radial gas distribu- 
tion. 

There are a few assumptions inherent to the model 
that need further justification: 

The assumed cooling function. The cooling function below 
10^ K depends on the degree of ionization of the gas. We 
made a few test runs with the cooling function of Spaans 
& Norman H1997|l truncated at the lower temperatures. 
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i.e. the cooling was set to zero at T < 10^ K, and found 
that the influence of the cooling at T < 10"* K on our re- 
sults was minimal, which is due to the fact that thermal 
energy of the shell is mostly lost at higher temperatures 
as mentioned earlier. 

Injection of energy by supernova explosions. In the present 
simulations we use a constant wind approximation (Mac 
Low & Ferrara I1999|l . However, our simulations have 
shown that the injection of thermal energy due to dis- 
crete SN explosions produce similar results. Indeed, the 
frequency of SN explosions in an instantaneously born 
stellar cluster remains nearly constant, if a Salpeter ini- 
tial mass function and an upper stellar mass of 100 Mq 
are assumed. This finding was also confirmed in numerical 
simulations of Mac Low & McCray l|1988|l . 
Star formation history. In the present simulations we have 
assumed a single starburst localized in the central region 
of Ho I that provides the energy for 30 Myr. Constant star 
formation in Ho I can be ruled out because of the low to- 
tal Ha luminosity of 4.3 x lO'^* ergs s~^ (Miller & Hodge 
I199()|l . which implies a SFR of only 0.004 Mq yr"^ Indeed, 
assuming the constant SFR of 0.004 AIq yr~^ and an age 
of Ho I of 10 Gyr, one obtains an estimate on the total 
stellar mass of 4 x 1O''M0, which is 2.5 times lower than 
that actually measured for Ho I (Ott et al. l2001|l . Variable 
star formation is plausible. However, in dirr's variations 
have a burst-hke nature (see e.g. Searle et al. 11973(1 . Our 
model corresponds to an explosive energy input following 
evolution of a single OB association with a total mass of 
3 X 10^ Mq . A comparison of the photometry presented 
in Ott et al. pMI with the STARBURST99 synthesis 
models (Leitherer et al. ll999|l shows that the B-band mag- 
nitude as well as the U— B color of the optical emission 
within the giant HI hole is in agreement with a stellar 
cluster 55 Myr of age (predicted by model 1) and a mass 
of ~ 2 X IO^Mq (Salpeter IMF, Z=0.004 metaUicity, upper 
and lower mass cutoff: 100 Mq and 1 Mq, respectively). 
Such a stellar population can easily provide the mechan- 
ical energy input to Ho I needed to satisfy model 1. The 
discrepancy of the observed stellar mass with our simula- 
tions (about 1 order of magnitude) is most likely due to 
confusion with older stellar populations. 

Nevertheless, we examine if variable star formation can 
affect our conclusions. We assume that the mechanical lu- 
minosity produced by SN explosions is time-dependent, 
for instance with a Gaussian shape. Its maximum is cen- 
tered at 15 Myr after the beginning of the energy input 
phase, with a FWHM of 8.5 Myr. This time-dependent 
mechanical luminosity is supposed to mimic the energy 
input from successively born stellar clusters within a sin- 
gle stellar association. To be consistent with the simula- 
tions in Sect. 01 we set the energy input phase to last for 
30 Myr and the total energy input equivalent to 300 and 
150 SNe for model 1 and model 2, respectively. Figure ITTI 
shows the azimuthally averaged gas distributions obtained 
at the same phases of the shell expansion as in Figs. |H1 and 
1101 As is obviously seen, model 1 with variable star forma- 
tion reproduces the observed distribution and the contrast 
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Fig. 11. Azimuthally averaged radial gas distribution ob- 
tained for a time-varying mechanical luminosity with a) 
model 1 and b) model 2 at three different phases of the 
shell expansion. The measured radial gas distribution of 
Ho I is plotted with the filled squares. 

in Sg between the central depression and the ring (the 
solid line, Fig. lllh .) even better than the same model with 
a single starburst approximation (the solid line, Fig. ISJ. 
Model 2 clearly fails to reproduce the measured contrast 
in Eg for both star formation frameworks. Our simulations 
also show that multiple, successively born smaller clusters 
are less effective in creating a bubble structure than a sin- 
gle starburst with the same total ejected energy. 
Numerical resolution. We made a few test runs with a 
higher resolution of 5 pc. Although the higher resolution 
simulations provide more details on the gas flow dynam- 
ics and the development of Rayleigh- Taylor and Kelvin- 
Helmholtz instabilities (Figs 01 and O, the integrated im- 
ages such as the HI spectrum, the channel maps, and the 
radial gas distributions (Figs.|HIEI and lllfl remain very 
similar to those derived with a lower resolution of 10 pc. 

6. Summary 

Numerical hydrodynamical modelling of the SN-driven 
shell formation is performed with the aim to reproduce 
the ring-like morphology of HI in the dirr Ho I. Two mod- 
els of Ho I are considered: a fast rotating galaxy with a 
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massive halo of 6.0 x 10^ Mq (model 1), and a slowly rotat- 
ing galaxy with a small halo of 4.0 x 10* Mq (model 2). 
Simulations reveal the following differences between the 
two models: 

1. The mechanical energy to create the ring of 1.7 kpc 

diameter in model 1, (3.0 ± 0.5) x lO''^^ ^j-gg qj. 399 ± 50 
Type II SNe, is a factor of 2 higher than in model 2. This is 
mainly due to the smaller gas volume density and weaker 
gravitational potential in model 2 as compared to model 1. 

2. In model 1 the shell breaks out of the disk before it 
expands out to its present size of 1.7 kpc. On the contrary, 
in model 2 the shell never breaks out of the disk, mostly 
due to the thicker vertical gas distribution and lower en- 
ergy input. As a consequence, both the measured contrast 
(^ 15) in gas surface density between the central depres- 
sion as well as the ring and the lack of gas expansion in 
the center of the ring are better reproduced in a model 
with a massive (6.0 x 10^ Mq) dark matter halo. 

3. Within model 1 an upper (i = 20°) and lower (i = 
15°) limit to the inclination of Ho I can be set from the 
modelled HI spectrum and the channel maps, respectively. 
This is probably the only handle to assess the inclination 
of gas-rich, nearly face-on dwarf galaxies. 

Our numerical hydrodynamical simulations suggest 
that Ho I was initially a halo-dominated dirr with an expo- 
nentially declining radial gas distribution. Approximately 
60 Myr ago a stellar cluster was born near the dynam- 
ical center of the galaxy from a parent gas cloud of 
(1.5 — 3.0) X 10^ Mq. Subsequent SN explosions created 
an expanding shell filled with the hot SN ejecta. About 
15 Myr ago the shell broke out of the disk pumping the 
hot gas into the halo. The fraction of the total gas mass of 
Ho I lost in a blow-out is low, within 1 — 2%. At present 
the shell has attained its biggest size of 1.7 kpc. In an- 
other 20 Myr the central hole would fill in and Ho I would 
appear as a dIrr, with a wave-like declining radial gas dis- 
tribution. 
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